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Abstract 

We use gauge/gravity duality to study deeply virtual Compton scattering (DVCS) 
in the limit of high center of mass energy at fixed momentum transfer, corresponding 
to the limit of low Bjorken x, where the process is dominated by the exchange of 
the pomeron. Using conformal Regge theory we review the form of the amplitude for 
pomeron exchange, both at strong and weak 't Hooft coupling. At strong coupling, 
the pomeron is described as the graviton Regge trajectory in AdS space, with a hard 
wall to mimic confinement effects. This model agrees with HERA data in a large 
kinematical range. The behavior of the DVCS cross section for very high energies, 
inside saturation, can be explained by a simple AdS black disk model. In a restricted 
kinematical window, this model agrees with HERA data as well. 
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1 Introduction 



The approximate conformal symmetry of QCD at high energies gives an ideal ground to test 
gauge/gravity duality [TJ in situations where there is available data from particle accelera- 
tors. Ideally one should be able to reconstruct the dual of QCD starting from its conformal 
limit, where the dual geometry is simply AdS space. The fact that the high energy limit of 
strings in AdS space does reproduce the observed hard scattering behaviour of QCD [2], as 
opposed to that of strings in flat space, is a clear indication that string theory on AdS is 
indeed a theory of the strong interactions. In this research programme, deep inelastic scat- 
tering (DIS) in QCD is among the processes that have been extensively studied using the 
AdS/CFT duality [3-28]. In particular, there has been progress in reproducing experimental 
data from HERA in the regime of low Bjorken x, where the gluons dominate the partonic 
distribution of the target proton [7][23]- This regime corresponds to the high energy Regge 
limit where the interaction between the probe and the target is dominated by pomeron 
exchange |29j . 

In this paper we shall test gauge/gravity duality models for pomeron exchange in QCD 
against the process of deeply virtual Compton scattering (DVCS), where an electron emits 
a space-like off-shell photon that scatters off a scalar target, and is then converted to an 
on-shell photon. Figure [T] represents this basic process. We will be able to match data from 
ZEUS [30] and HI [31] experiments with very good results. Ultimately we want to make 
these models a useful tool to match experiments and to give new predictions for QCD. There 
have been some works that considered DVCS in gauge/gravity duality [32], [33], [33] - Here we 
make direct link to experiment, following our earlier work [?J 123] ■ 

DVCS has a long history of study from QCD. Although a review of this approach is 
outside the scope of this paper (see for example [33] and references therein, or [SB] for a 
quick introduction), let us briefly comment on some recent results that relate to pomeron 
exchange. The authors of [37] use the dipole model to fit the HERA data for exclusive 
diffractive processes. In this picture the incoming off-shell photon fluctuates into a quark- 
antiquark dipole, which then interacts with the proton via the exchange of a gluon ladder. 
They do not use DVCS data to fit their parameters, but still get good agreement when 
they compare it to their model. A recent paper [3H] fits the DVCS (and vector meson 
production) data using a non-linear pomeron trajectory, with good results. However, we 
feel an advantage of our approach is that due to our strong coupling starting point we have 
better access to the non-perturbative region, and due to the diffusion in AdS space of our 
pomeron kernel we can fit the data for all of the Q 2 values, i.e. in both the soft and the hard 
scattering regions, using just a single set of parameters. Furthermore, we fit the differential 
cross section directly as a function of kinematic invariants Q 2 , W and t, and the cross section 
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Figure 1: In DVCS the electron couples with a basic QED vertex. We focus on the QCD process 
7*p —7- 7p, dominated at weak coupling by the exchange of the BFKL pomeron. There are three 
kinematical invariants: the offshellness of the incoming photon Q 2 ; the center of mass energy 
W = y/s; the momentum transfer 

as a function of Q 2 and W, without the need to fix one of the variables and only fit the 
dependence on the other one. 

In section [2] we revisit in some detail gauge/gravity models for pomeron exchange in 
QCD. Then in sections [3] and [4] we do the computations of the cross section for 7*p — > jp, 
and compare with data in section [5j Section [6] concludes. 

2 Review of pomeron in gauge/gravity duality 

In many 2—7-2 QCD processes the Regge behaviour 



where s and t are the Mandelstam variables, with s 3> t, is observed. The exchanged state 
can be thought of as a reggeon, with the intercept jo — a (0) determining the growth of total 
cross sections, a to t ~ s j0_1 . In particular, this behaviour is observed in processes for which 
a pomeron - the leading state with the quantum numbers of the vacuum - is exchanged. 
In perturbative QCD, this story can be made very explicit by considering the hard BFKL 
pomeron, which is a perturbative resummation in (A Ins)™, where A is the 't Hooft coupling 
The hard pomeron describes the exchange of two reggeized gluons, with many ladder 
diagrams also made of reggeized gluons, and contains only planar diagrams, so that it is 
the dominant contribution in a large N expansion (see for instance figure [If. For processes 
where the QCD scale is important, for instance whenever we probe a target proton at the 
QCD scale, Regge behaviour is still observed, though with a different value of the intercept 



A(s,t) ~ 8' 



a(t) 




(1) 
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jo- In this case the exchanged state is usually called the soft pomeron and it dominates the 
total cross section. 

Gauge/gravity duality gives a beautiful new insight into the above discussion |4j. The 
main idea is that the hard and soft pomerons correspond actually to the very same Regge 
trajectory, with intercept depending on the effective coupling constant, based on the AdS 
radial coordinate. The pomeron Regge trajectory is identified with the BFKL hard pomeron 
at weak coupling, while at strong coupling it is identified with the reggeized graviton tra- 
jectory of strings in AdS. As the coupling is increased, the intercept j = jo(\) varies from 
one to two. These ideas can be made very precise in M = 4 SYM, and suggest that one 
can analyse processes dominated by pomeron exchange starting from the conformal limit of 
QCD j7]. In this limit, we have the BFKL hard pomeron as the weak coupling description. 
However, now we also have a strong coupling conformal description, from which we can 
perturb around. More precisely, we may consider a reggeized graviton in AdS§ space and 
analyse its predictions for experimental data in regions where strong coupling effects are 
important [24] . 

Our starting point is QCD at some fixed value of the coupling (in DIS and DVCS, 
for instance, this is defined by the scale of the probe given by the off shell photon), in a 
conformal kinematical window where all the scales are above the QCD scale. Let us remark 
that in this window perturbation theory does not necessarily work, since there may be a 
kinematical enhancement of the effective coupling. This is in fact the case of low x physics, 
where non-linear effects give rise to a large gluonic component in the parton distribution 
functions. 

In the conformal limit the momentum space amplitude A(ki) can be expressed in terms 
of a reduced conformal amplitude B(S, L), which depends only on two conformal cross ratios 
5* and L, that have a natural interpretation in terms of a dual AdS scattering process. To 
see how this arises start from the four point function of the electromagnetic current operator 
j a {y) of dimension £ = 3 and of a scalar operator 0(y) of dimension A, 



For simplicity we suppress, in this section, the indices of the current operator. The full 
computation is presented in detail in [17J. The Regge limit corresponds to the limit where 
the points are sent to null infinity. More concretely, defining light-cone coordinates 
V = (y + iV - )2/J-)) where y± is a point in transverse space M 2 , this limit is attained by 
sending 

y\ ->■ -oo , yt -> +oo , y 2 -> -oo , yl -> +oo , 
while keeping yf and yi± fixed. This limit is better studied if we do conformal transforma- 




(2) 
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tions on the points y^ according to 

1 



x i =(xf,x i ,x iL ) = — x (l,y?,2/ii) , * = 1,3, (3) 

Vi 

Xi = (xf, xj, x i± ) = — - (1, y" 2 , y i± ) , i = 2, 4 . (4) 

Vi 

The conformally transformed amplitude 

A(x t ) = (j(x 1 )0(x 2 )j(x 3 )0(x 4 )) , (5) 

is related to (|2j) by the standard conformal transformation rules for primary operators. One 
can then see that in the Regge limit the amplitude depends only on the small vectors 

x = X\ — x 3 , x = x 2 — x 4 • (6) 

Moreover, since conformal symmetry imposes that the amplitude is a function of two cross 
ratios, we can write the only two independent invariants as 

a 2 = x 2 x 2 , cosh p = — X X , (7) 

\x\\x\ 

so that the Regge limit corresponds to sending cr — ^ with p fixed. Finally, we can introduce 
the Fourier transform [33] 

A(x,x) = J dpdpe- 2ip - x - 2m B{p,p) . (8) 

A careful analysis of the ie prescription for the amplitude A(x,x) shows that B(p,p) only 
has support inside the future light cone. We may then take p and p as future directed 
timelike vectors, and show that B(p,p) has the form 

B(S,L) 



B(p,p) 



(-p 2 ) 2 -5(-p2)2-A 



The reduced amplitude B(S, L) is extensively discussed bellow, and depends only on the 
two cross ratios 

S = A\p\\p\ } coshL = -^|. (9) 

The Regge limit corresponds now to large S and fixed L. 

The sequence of transformations, just described, are schematically represented in figure 
[2j The equation relating A(ki) to B(S, L) will only be needed in section [3j for now we want 
to discuss the form of the amplitude B(S, L) when it is dominated by the exchange of a 
reggeon, as dictated by conformal Regge theory [40J. That the Regge limit of large s and 



5 



FT. 

Figure 2: As usual the amplitudes A{ki) and A(yi) are related by Fourier transform (F.T.). In the 
Regge limit we consider conformal transformations (C.T.) of the points yi, defining an amplitude 
A(x,x). With a further Fourier transformation we introduce the amplitude B(p,p). The relation 
between the original amplitude A(k{) and B(p,p) is the one we would write using AdS/CFT, with 
the specific function for B(p,p) depending on the corresponding Witten diagram. 

fixed t corresponds to the kinematic limit of large S and fixed L, will also become clear when 
we directly relate the cross ratios S and L to the Maldelstam variable s and the impact 
parameter l±, in section [3] Finally, note that the cross ratios have a natural interpretation 
from the view point of a dual AdS scattering process, where the transverse space is a three- 
dimensional hyperbolic space if 3, whose boundary is conformal to the physical transverse 
space IR 2 . The cross ratio L is then identified with the geodesic distance between two points 
in if 3, and can be related to the physical impact parameter l± along IR 2 . The other cross 
ratio S measures the local energy squared of the scattering process in AdS, which justifies 
its role as the large energy in the conformal Regge theory reviewed below. 

In order to associate a scattering matrix involving a target hadron to the conformal 
amplitude B(S, L), we need to introduce a discrete spectrum in the theory therefore breaking 
conformal symmetry. AdS/CFT will instruct us how to do this, without changing the 
conformal structure of B(S, L). This will be the case of the AdS black disk and conformal 
pomeron models considered in this paper. To include the effects of confinement in the 
amplitude B(S, L), we will then consider the hard wall pomeron model. For another recent 
work that explores the conformal symmetry in the context of low x DIS see |41j . 

2.1 Regge Theory in CFTs 

Let us now see how conformal Regge theory [ID] restricts the form of the reduced conformal 
amplitude B(S, L) introduced above. We shall assume that the amplitude for low x DVCS 
is dominated by a Regge pole, corresponding to the QCD pomeron. This Regge pole is 



C.T. 

< »• 



AdS/CFT 



FT. 
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associated to the exchange of the twist two operators in the leading Regge trajectory. The 
contribution of all these exchanges can then be written as 

B(S, L) = ^J2( SJ + C" 5 ) ') / du h ^ n ^ ■ ( 10 ) 

J>2 J 

This form of the amplitude arises from the conformal partial wave expansion of the ampli- 
tude. In particular, for fixed spin J of the exchanged operator, the S dependence is given by 
if?" 7-1 , and the dependence in the impact parameter L is determined by the coupling of the 
external operators to the exchanged operator and also by the dimension of this exchanged 
operator, which appears as a pole of hj(u) at iu = A (J) — 2. The function 

^> = it^ <»> 

is a basis of harmonic functions on a three-dimensional hyperbolic space if 3. For two points 
on x, x G H 3 , L is defined as the geodesic distance between these two points. A simple 
calculation shows that the function Q iu satisfies (O x + 1 + u 2 ) Q iu (x,x) = 0, and that 

J duVt iu (x,x) = 5 Hz (x,x) . (12) 



To compute the sum in (10) one re-writes it as a Sommerfeld- Watson integral over spin 

J, 

7T r dj s J + {-sy 



B(S ' L > = -2S J 2^^7T J ^WW), (13) 
where the integration contour in J arises from picking all the contributions from the twist 
two operators of spins J = 2,4, • • • . Deforming the J-contour one picks a Regge pole at 
J = j(v), which gives the leading behaviour of the amplitude in the Regge limit 

B(S, L) = — J dufiiy) S^- 1 nUL) , (14) 

Both functions j(y) and fi{v) are also functions of the coupling constant. The Regge spin 
j{y) is related, by the analytic continuation A = 2 + iu, to the spin-anomalous dimension 
curve A = A(j) of the twist two operators on the leading Regge trajectory. From the weak 
coupling side, both in QCD and in TV" = 4 SYM, we know that 



m x) = 1 + A _ * (l±i^) _ * (i^)) + oM , (15) 

where A = ^y M A^ is the 't Hooft coupling. On the other hand, from general arguments, we 
expect that string theory in AdS§ will give the following strong coupling expansion for a 
general class of conformal theories jl] 

j(u,\) = 2-±±£ + 0{l/\), (16) 
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where here the 't Hooft coupling A = R A /ot' 2 is defined by the ratio of AdS radius to the 
dual string a', and obeys A ~ Qym^i although the overall numerical factor is only know in 
specific examples such as TV = 4 SYM. 

It is straightforward to relate the function fl{p) to the residue of hj(v) in the complex 
J-plane. The function /3(f) is determined by the reggeon exchange and by the coupling of 
the reggeon to the external states, and can be written in the form 

$(v) = W(v)Q(v)W(v). (17) 

The functions W(f) and W(f) are impact factors for the external operators [12], and were 
computed at weak coupling for the QCD current operator j M = ip^^ip in [17]. At weak 
coupling the impact factors can also be directly related to the dipole wave functions of the 
probe and target [7J. In M = 4 SYM the next-to-leading order corrections of the impact 
factors were computed in [43]. These functions can also be computed at strong coupling by 
assuming minimal coupling between external states and the exchanged string states in AdS 
[40J. The function Q(f) is related to the exchanged reggeon. 

To better understand the Regge residue it is convenient to write the amplitude B(S, L) 
in the original form, as a sum over the spin of the exchanged twist two operators, 

^ L >-^I^ s Wr B(J ' Ll (18) 

Equation (13) gives the radial Fourier decomposition, 

B(J,L) = J dfhj(f)n iv (L). (19) 
Near the Regge pole J ~ j(f) the Fourier component hj(f) will be of the form 

27T? 

hj(f) = I(u) G(J, f) 1(f) . (20) 



The impact factors 1(f) and 1(f) can be simply related to those in (17), after computing 



the J integral. The dependence on J of G( J, v) is given precisely by the Regge pole 

G(J, v) = , . (21) 

u*+(A(J)-2) 2 

For the purpose of computing the J integral by picking the pole at J = j(f), we only need 
to determine the leading behaviour of this function near the Regge pole, which is defined 
by the condition ±iu = A(j(u)) — 2. Thus, near the pole we can write 

G( J ' U) = 9- M-( l(l '( W ' (22) 

2tfA'(j(f))(J-j(f)) 
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where A' (J) is the derivative of the anomalous dimension function A = A (J) for the twist 
two operators in the leading Regge trajectory. At strong coupling this function reduces to 
the familiar form [I] 

G ^»mhw)' (23) 

but ( 22 ) is valid for any value of the coupling. 

The function G(J, v) that gives rise to the Regge pole can be related to the reggeon 
propagator in AdS space. To understand this we consider again the radial Fourier decom- 
position 

G(J,L) = J dvG(J,v)tt iu (L) . (24) 

This function is a propagator in the transverse space H 3 of the dual AdS scattering process, 
satisfying the propagator equation 

n H3 - 3 - A (J) (A(J) - 4)] G(J, L) = -5 H3 (x, x) , (25) 

where we recall that L is the geodesic distance between x and x in H$. 

Finally we wish to understand how the above H3 propagator can arrive as the exchange 
in AdS of a field of spin J and dimension A. When computing the Witten diagram for the 
exchange of such a field between two external particles in the Regge limit, one learns that 
the amplitude depends on the integrated propagator of the exchanged field along the null 
geodesies of each external particles. More precisely, from the propagator of the exchanged 
field one obtains 

/oo 
dXdX n^(X(A),X(A)) = -2iS J - 1 G{J ) L) , (26) 
■00 

where here A and A are the affine parameters along the particles geodesies X = X(X) and 
X = X(X). In this equation the invariant S is given by S = —2K ■ K, with K and K the 
tangent vectors to the geodesies. The propagator lT J ) is the contraction of these tangent 
vectors with the spin J propagator, i.e. 

n( J )(X,X) = (-2) J K^---K aj K^---K^ 1I,„... ( , / ;,.. ;; (.Y..Y). (27) 



One can then show [44J that equation (25) for the H3 propagator arises from the propagator 
equation for a spin J field of dimension A that, in the Regge limit, is given by 

n AdSs - A(J)(A(J) - 4) + J n ai ... aj p v .. Pj (X,X) = igaxifc ■ ■■9\aj\l3 J )SAdS 5 (X,X) + ■■■ , 

(28) 

where ■ ■ ■ represent terms that are negligible. For instance, in the infinite coupling limit, 
where only the graviton contributes to the exchange, by setting A(2) = 4 one obtains pre- 
cisely the propagator that appears in the Regge limit of the corresponding Witten diagram 
for graviton exchange 
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2.2 Pomeron at strong coupling 



Let us finally review the computation of the amplitude B(S, L) for pomeron exchange at 
strong coupling |U 46] Q First we compute the v integral in (19) and then the J integral in 
(18). In particular we are interested in the strong coupling expansion of the reggeon spin 
given by (16). Using the explicit expressions for the harmonic functions Qi U in (11) and for 
the pomeron propagator in (21), the v integral can be done by residues with the result 

-ivh 



(29) 



where in this expression v = v{J) = —i(A(J) — 2). To perform the J integral we need to 
consider first the spin-anomalous dimension relation A = A( J), with inverse J = j(y) given 
in (15) or (16). We shall see that the J integral can be done by expanding in v. Both at 
weak and strong coupling we can write 

j{v)=j Q -Vv 2 + -- - , (30) 

In particular, at strong coupling and keeping only the terms in 1/a/A, the v expansion stops 
precisely at order v 2 [4"0~] . with 



Jo = 2 - 



V 



2VA 



(31) 



Thus, in general the v dependence in (29) is given by 



iv(J) = A(J) - 2 = -i 



(jo ~ J) 
V 



(32) 



Hence, to perform the J integral, we note that the integrand (29) will have a branch cut for 
Re (J) < jo, arising from the argument of the exponential. Indeed the impact factors are 
even functions of z/, admitting an expansion in powers of u 2 , so they do not give rise to any 
additional branch cut. After a simple analysis of the contour of integration we can re-write 
dll as 



B(S,L) 



N 2 4vrsinhL S 



l0 dJ SJ + 



-sy 



sm(vL) I{u) I{v) . 



(33) 



sin(7rJ) 

Next, using the relation v = v{J) in (32), we change from the integration variable J to the 
variable z/, with the result 



B(S, L) = - 



1 vs jo ~ 



N 2 47rsinhL 



du v 1 + 6 e -(\nS)V^-iuL j,. J, x 

sin(7rj(i/)) 



(34) 



^^Regge trajectories at strong coupling of other states, for example the odderon, can be similarly found 
(see [57] for details). 
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where we note that j = j(y) is given by (30) 



We are interested in the high energy limit, at strong coupling, with the exchange of all 
twist two operators in the leading Regge trajectory. This means that we should consider 
the limit of Pin 5* = \a.S/{2y/X) ^> 1 (this is not the strict A — > oo limit, where only the 
energy momentum tensor contributes to the exchange). The previous integral can be done 
using the saddle point approximation]^] The saddle point is at 

Uo = 2^S ' (35) 

which is close to zero, so it is consistent with an expansion in v. A straightforward compu- 
tation gives 



where In = 1(0) is the intercept value for the impact factor I(y) and similarly for I(v). In 
addition to the expansion in a/X/ In S <C 1 used to compute the J integral, we are also doing 
a large 't Hooft coupling expansion in 1/VA <C 1. 

To close this section we consider the overall normalisation of the amplitude, in particular 
the dependence on the 't Hooft coupling. So far we took S as a dimensionless cross ratio. 
When we consider a specific conformal field theory with an AdS string dual there are two 
new length scales, the AdS radius R, and the string tension as defined by a'. The AdS 
radius allows us to define a dimensionful S, whose physical interpretation is related to the 
energy squared with respect to global AdS time. This scale is enough to reproduce the 
amplitude given by graviton exchange in AdS from the above formulae, and no a' enters 
the final result. In this case the amplitude can be obtained simply by doing the integral in 



(34) by picking the poles at the zeros of sin(7rj(i/)) — 7r(2 — jo) — nT>v 2 . The result is 

2L 



i ^ 2r e 



BiS^^-miMR'S^—^. (37) 

Here S has dimension of energy squared, as determined by the AdS radius R, so we replaced 
in this expression the original dimensionless S in (34) by R 2 S, since the conformal amplitude 



B(S, L) is dimensionless. Looking at the gravity computation done in appendix [A] where 
both for external currents and scalar fields we have 



2 Note that in the weak coupling limit, we must require Din S « A In S 1 to compute this integral 
using the same saddle point approximation. 
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we can read off the impact factors, and in particular obtain the strong coupling result 
J(2i) = J(2i) = 7T. 

On the other hand, if we consider the graviton Regge trajectory derived from strings in 
AdS as in (36), the string tension plays an important role. Specifically, after introducing 
the AdS length so that S is dimensionful, the dimensionless combination that enters in 
the above expressions in the place of S, say at the very beginning in equation (18) for the 
amplitude, is precisely a'S. Implementing this modification, and normalizing the overall 
amplitude such that it gives the correct gravity result in the limit A — > oo, as computed 
above, equation (18) is correctly normalized by 

dJ (a'S) J H 



B(S,L) 



K 2 



71 



a' 2(a'S) 

This means that we should replace S by 



-a'sy 



a'S 



lixi sin(7rJ) 

zzs 



B(J,L). 



(39) 



(40) 



in the expression (36), and also multiply this equation by an overall factor R 2 /a' . Thus, 
our final result for the strong coupling gravi-reggeon amplitude is 



Tip 

Y 



B(S,L) = gl (1 + 2 cot 
where the parameters p and g$ are defined by 

p = 2 - j = AV 



a'S) 



pin(a'S) 



(pln(a^)) 3 / 2 sinhL ' 



2 
71 



9l 



1 hh 



(41) 



(42) 



2^ A 2 

and we considered expansions in 1/y/X <C 1 and vA/lnS 1 <C 1. The derivation in this 
section is similar to the one in [46 J , but the way we are presenting the result directly relates 
the coupling g\ to the impact factors of the external states. 

Let us finally comment on the weak coupling limit of the amplitude, corresponding to 
the exchange of the hard BFKL pomeron. The very same computation leads again to the 
result (36), just that this time jo ~ 1 + A%(0) and V w — Ax"(0)/2, where the function x{ v ) 
can be read off from (15). Taking care of the overall normalization, as in (39), we obtain 

{a'sy^ 1 



B(S,L) 



So 2 



XV 1 



-(?) 



g 4Bln(a'S) 



{AV\n{a'S)fl 2 sinhL 



(43) 



where j is the weak coupling intercept, g 2 is still related to the impact factors by (42) and 
(a'S) is given by (40). This expression is obtained by considering the expansions in A 1 
and l/(AlnS) 1. For instance, we can consider the amplitude between two dipoles. In 
this case each impact factor comes multiplied by A, as can be seen from the diagram of 
figure [l] in the introduction, and as computed explicitly in [42j [17] . This computation gives 
therefore the correct leading term at order A 2 , which is real. 
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Figure 3: Conformal amplitude B(S, L) for fixed S obtained from exponentiating the single reggeon 



exchange. The picture was obtained from the expression (41 ) for the reggeized graviton trajectory. 
By varying S it is also simple to verify that the size of the black disk grows linearly with In S. 



2.3 AdS black disk from pomeron exchange 

Both at weak and strong coupling, one can study the dependence of the amplitude B(S, L) 
with L, for fixed S. In both cases the amplitude reaches unity for some value of L = L S (S). 
One then needs to include corrections that come from both multi-ladder reggeon exchange 
and also from reggeon-reggeon couplings. We may then define the amplitude in terms of a 
phase shift x(S, L) such that 

B(S, L) = 1 - e ix{s ' L) . (44) 

Provided that the imaginary part of the phase shift grows fast enough as L decreases, the full 
amplitude will in general have the form given in figure |3j This suggests the introduction of a 
conformal (AdS) black disk model for the amplitude, as also plotted in figure [3] The growth 
of the disk with energy can be determined by defining first Imx(>S', L s ) ~ 1. Approximating 
the full amplitude by the single reggeon exchange ( 14 ) , one can then compute the v integral 
by saddle point, to obtain the general behaviour L s = uhxS, for some constant u. This 
computation can be done very explicitly both at weak coupling [7] , for the BFKL pomeron, 
or at strong coupling, for the BPST pomeron (reggeized AdS graviton). 

For the AdS black disk model we simply assume that the imaginary part of the phase 
shift grows fast enough with decreasing impact parameter L, therefore leading to a black 
disk, independently of the details of the phase shift and of the value of the coupling constant. 
For instance, figure [3] uses for the phase shift the amplitude for the exchange of a single 
AdS pomeron ( |4l| ). Simple combinatorics show that exponentiating the single exchange 
amplitude corresponds to the inclusion of all ladder reggeon exchanges, while dropping 
reggeon couplings 
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2.4 Conformal symmetry breaking 

This paper considers exclusive photon production in the Regge limit, when a space-like 
off-shell photon scatters off a scalar target, with the final state given by the target and an 
on-shell photon, as represented in figure [T] in the Introduction. As usual the Regge limit is 
defined by taking large s = — (£a + A>2) 2 with fixed momentum transfer t = — (ki + k^) 2 = —q\ 
and virtualities k 2 . 

As already explained we shall be interested in the kinematical regime where QCD is 
approximately conformal. At weak coupling this corresponds to the limit where the exchange 
is dominated by the BFKL hard pomeron, while at strong coupling it corresponds to the 
exchange of the graviton Regge trajectory in the dual AdS space. To be more precise, our 
starting point is the conformal limit of QCD. We then break conformal symmetry using the 
gauge/gravity duality by introducing a hard wall in the dual AdS space. Thus we start with 
the high energy Regge limit 

s > -t , Q 2 , (45) 

where Q 2 = k\ > is the off-shellness of the incoming photon. As usual the Bjorken variable 
is given by 

Q 2 ^Q 2 fAR . 
x = -2k~^ T' (46) 

where the last equality is valid in the high energy Regge limit of x <C 1 here considered. 

Conformal invariance requires that all the above energy invariants are above the QCD 
scale Aqcd ~ 0.2 GeV. In particular, the condition Q ^> Aq CD means that the photon 
probe is much smaller that the size of the proton, so that corrections in Aqcd/Q that break 
conformal invariance are suppressed. Most data used in this paper satisfies this condition 
(some data has Q ~ 1.5 GeV, which is still a few times larger than Aqcd)- Moreover, we 
must also satisfy the condition —t ^> A 2 q CD for the momentum transfer, so that corrections 
in A 2 Q CD /t are also suppressed. Again this condition is satisfied for most of the data here 
considered {—t will range from 0.1 to 1 GeV 2 , so that we will be above A 2 q CD ~ 0.04 GeV 2 , 
and only for the lowest value of —t will the approximate conformal symmetry begin to show 
signs of breaking down). 

Another source for conformal symmetry breaking comes from the masses of quarks, which 
also introduce new scales. To see this note that the coupling between the off-shell photon 
and the target is determined by the electromagnetic current operator j a = e$^/ Qf^fja^f, 
where the sum is over quark flavours. The off-shell photon creates a quark/anti-quark pair 
which then couples to the target through the exchange of many gluons. We therefore need 
to work in a range of Q such that 2m f <C Q. For these flavours the quarks contribute to 
the process and can be taken as massless. For the other flavours we need to be below the 



14 



scale defined by the threshold for pair creation, i.e. Q <C 2m./, so that these quark loops 
can be ignored. We shall be working in a kinematical region for which 

1.5 £Q£ 8.4 GeV. (47) 

Therefore we can safely take the u, d and s quarks as massless. For the charm quark, with 
mass m c = 1.27 GeV, some of the data points will be in the region where its mass can 
not be neglected. Since the majority of data points have Q above the 2m q threshold we 
will assume this quark enters the process and is massless (for the analysis of differential 
cross section there will be 4 points below the threshold out of 52 , and for the total cross 
section 13 points out of 44). This is clearly an approximation that enters the assumption 
of conformal symmetry^ Finally, for the bottom quark, we can ignore its existence, except 
for a single data point at Q — 8.4 GeV, since m& = 4.2 GeV. Thus, we shall analyse a 
kinematical region where the four quarks u, d, s and c are approximately massless and the 
remaining ones can be neglected. 

Both the AdS black disk and conformal pomeron models considered in this paper, assume 
that conformal symmetry dictates the form of the amplitude. The very good matching of 
experimental data obtained for both models clearly points out that we are indeed analysing 
a conformal window of QCD. Conformal symmetry breaking only enters through the as- 
sumption that the theory contains a discrete spectrum of states, as necessary to define a 
target hadron. For the black disk, the final expression is independent of the coupling of 
the reggeon to the external states, so that the model depends only on the parameter u 



introduced in section |2.3[ on a scale defined by the target hadron, and on a overall con- 
stant determined by the normalization of the current operator two-point function. This 
model will match data very successfully in a kinematical window where so far there are 
only a few available experimental points. For the conformal pomeron, we need to include 
the coupling of the reggeon to external states. This is one of the parameters of this model 
that multiplies the normalization of the current operator two-point function. Together with 
the scale of the target hadron and the value of the intercept jo, this model also has three 
independent parameters. We shall see that the conformal pomeron matches data in a quite 
broad kinematical range. 

To include QCD effects associated to confinement, we proceed phenomeno logically by 
considering the AdS hard wall pomeron model, where the dual AdS space finishes at some 
fixed scale, defining an additional fitting parameter, which should be around Aq C d- In 
this phenomenological model one is then able to take into account the dependence of the 
amplitude on new dimensionless ratios constructed from Aqcd, like Q/Aqcd and tj 'Aq CD , 



3 It also enters the weak coupling computation of the normalization of the currrent operator explained 



in section 3.2.1 but this plays no role in the fit to experimental data. 



15 



which break conformal invariance (as it is usually the case in gauge/gravity models of QCD, 
this is a phenomenological model, since we do not have control in terms of a well defined 
expansion in these small parameters). Thus we expect this model to be more accurate in 
regions where either the scale set by the probe, or the momentum transfer, become close to 
h-QCD- The fit to DVCS data will confirm these expectations. 



3 Hadronic Compton tensor in the Regge limit 

In this section we shall relate the results presented above for the conformal amplitude 
B(S,L), given by the conformal Regge theory, to the amplitude for the 2—^2 process in 
DVCS, as represented in figure [T] This uses the conformal impact parameter representation 
introduced in [2H] and extensively discussed in [T7] . We then compute the hadronic Compton 
tensor, which is the basic object that allows one to compute the 7*p — > -yp cross section. 

Before starting let us note that in the explicit computations here presented we use light- 
cone coordinates (+, — , _L), with metric given by ds 2 = —dx + dx~ + dx 2 ± , where x± e M 2 is 
a vector in the impact parameter space. For the incoming particles we take 

k^fyTs-^o] , k2= (Ml^ s ^o] t (48) 



where M is the mass of the target, and the incoming off-shell photon is space-like with 
k\ = Q 2 > 0. For the outgoing particles 

. ( r ll-Q' 2 \ h ( M 2 + ql \ 

fa = -[Vs, 7= — ,q± , «4 = - 7= — ,Vs, -q± )■ (49) 



For the outgoing photon we consider k\ = Q' 2 > 0, and then take the Q' — >■ on-shell limit. 



In section |3.2.4| we briefly consider lepton pair production for which the out-going photon 
is time-like. Note that the variable s, as defined in the above vectors, coincides with the 
Mandelstam s only in the Regge limit here considered. 



3.1 Conformal impact parameter representation 

As already discussed, we define the DVCS amplitude starting from a conformal field theory. 
Conformal symmetry breaking will then be implemented following the geometrical intuition 
gained from the AdS interpretation of the scattering amplitude. The basic object we need 
to consider is the momentum space correlation function 

(27r) 4 5 (J2 k)iT a \k 3 ) = (fihWk^fihWh)) , (50) 
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involving the electromagnetic vector current j a and a scalar primary O of dimension A. In 
the conformal Regge limit jlQ] , this correlation function can be written as [391 HZ] 

T a \k 3 ) « 2zs y dl x e iq ^ J J ^ ^(z) $ (^) t 1 - e ix(5 ' L) ] M T , (51) 

where the phase shift x^V * s a tensor that encodes all the dynamical information and, due 
to conformal symmetry, depends only on the variables 

^2 _|_ ^2 _|_ p 

S = zzs , cosh L = — . (52) 

2zz V ; 



In section 



2.1 



we refered to this conformal amplitude as B(S,L) = 1 — e %x ^ s,L \ and for 
simplicity we omitted the indices in that section. For a conserved current the Greek indices 
\i and t label tangent directions to a three-dimensional hyperbolic space if 3, with metric in 
Poincare coordinates given by 

m = + m 

We work in units such that the AdS radius R — 1, otherwise there would be an overall 
factor of R 2 in this metric, and the cross ratio S would be given instead by S = zzs/R 2 . 
The scalar function $ and the tensor function ^ ab J are associated with the operators O 
and j a . Their explicit form was given in |17j . 

The above conformal representation ( |51~| is quite natural from the view point of the dual 
AdS scattering process, where transverse space is precisely a three-dimensional hyperbolic 
space i^3, whose boundary is conformal to the physical transverse space M 2 . The cross ratio 
L is then identified with the geodesic distance between two points in H 3 that are separated 
by l± along IR 2 and have radial coordinates z and z. The other cross ratio S measures the 
local energy squared of the scattering process in AdS. The directions tangent to H3, which 
we label with Greek indices, are now the physical polarizations of the AdS gauge field dual 
to the vector operator j a . Moreover the functions ^^^(z) and $(2) can be related to the 
bulk to boundary propagators (non-normalizable modes produced by a plane wave source 
created by the dual operator at the boundary) of the dual AdS fields and their couplings 
to the exchanged reggeon, as shown in the strong coupling computation (i.e. AdS Witten 
diagram) presented in appendix [A} 

It is important to note that the above conformal representation of the amplitude is valid 
for any value of the coupling constant, since it relies only on conformal invariance. At weak 
coupling we can reproduce this exact form of the correlation function, with the phase shift 
determined by the hard pomeron exchange [17]. The functions ^ ah J(z) and &(z) can then 
be related to the dipole wave functions of the current and scalar operators [7J. At strong 
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coupling this form of the amplitude is reproduced by computing the AdS Witten diagram 
given by t-channel graviton exchange. This is a very non-trivial check of the all coupling 



form of the conformal impact representation (51). To verify this fact, and to convince the 



sceptical reader, we present this computation in appendix A| 



In the Regge limit the amplitude is dominated by the pomeron Regge trajectory. At 
leading order in M = 4 SYM, as well as in the A — > oo limit in A^ = 4 SYM, the phase shift 
X M r is diagonal. For simplicity, we assume that this holds for any value of the coupling and 
write, throughout this paper, 



■XT 



X(S, L) 8» T . 



(54) 



At strong coupling the phase shift is diagonal because the coupling between the gauge field 
and the graviton in AdS is determined by the energy-momentum tensor, which singles out 
the T ++ component in the Regge limit. Recalling the construction of the pomeron vertex 
operator [4] from the graviton Regge trajectory in AdS, we also expect the amplitude for 
this Regge trajectory to remain diagonal. 

Before we show how to compute the hadronic tensor that results from the impact param- 
eter representation, let us note that for a diagonal phase shift the contraction ty ab = ^f (ib ^{z) 



appearing in (51) has the explicit form 



/ (Q' 2 + q 2 ± )^ v^l^o \ 



Q 2 * lQHQ /2 + ql)Vo j^QW^t 



where 



V 



7T 



(55) 



QQ'S^t J 



y n = -C—z*K n (Qz)K n (Q'z), (n = 0,l) (56) 

with K n the modified Bessel function. For completeness let us also write the explicit form 



of the scalar function appearing in (51) 



r(A)r(A - 1) V 4 



QQ' 



A-2 



z 4 K 2 „ A (Qz)K 2 _ A (Q'z). 



(57) 



Although the reader may recognise in this expression the form of bulk to boundary AdS 
propagators, let us emphasise again that these expressions were obtained by simply assuming 
conformal symmetry, and are valid for any value of the coupling constant. Finally, note that 
the constants C and C are fixed by the normalization of the two-point functions for the 



4 The very same Witten diagram was computed in 
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external operators that give the disconnected piece of the four point function (50), with 

2 ab _ r } 1 A 1 b ft 



3.2 Hadronic tensor 

The hadronic Compton tensor is closely related to the above correlation function. It is given 
by the matrix element of the time-ordered product of two electromagnetic currents 

W a \k 3 ) = i J d*ye^(h\T {f(y)j b (0)} \h) , (59) 

where \ki) represents the hadronic state with momentum ki (i = 2,4). 



We can now use the impact parameter representation (51 ) to write the hadronic Compton 
tensor as 

/r A 7 A? 
dl ± e iq ^ I ^^* ab (z)<f>(z) [l- e lx{s ' L) ] , (60) 



where we already assumed a diagonal phase shift as explained in (54). The subtlety in this 



equation, in comparison with (51) for the correlation function (50), is in the function $(z) 
associated to the target. At strong coupling, as can be seen in the correlation function 
computation presented in appendix [A] this function is related to the product of two bulk 
to boundary propagators of the scalar field (p. The gauge/gravity duality gives then a clear 
geometric picture on how to break conformal symmetry, therefore introducing a discrete 
spectrum in the theory. This is done by deforming the AdS space in the IR region of large 
holographic variable z, and the simplest way to do this is to introduce a hard wall in the 
AdS space at z = z [2jj. Moreover, just as for strong coupling, the hard wall model can also 



be defined using the conformal impact parameter representation (51 ), valid for any value of 
the coupling. In particular, this may also define a way of breaking conformal symmetry for 
the weak coupling BFKL pomeron. We comment on this point in the Conclusion. 

We therefore introduce a IR cut off z , and follow our strong coupling AdS intuition. In 
this respect there are three aspects that need to be taken into account. Firstly, the scalar 
field has a discrete spectrum with normalizable modes <p n {z) satisfying 

dz 

— (f)* m (z)(f> n (z) = S m ,n- (61) 

The function *&(z) should therefore be determined by one such normalizable mode 4> n (z). 
In this paper we do not take into account the details of the target function, instead we shall 
use the simplest possible function, i.e. we consider the Dirac delta function 

$(z) = z 3 S{z - z m ) , (62) 
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where defines the scale of the target. This approximation will simplify considerably our 
formulas, still allowing for an excellent fit to experimental data. 

Secondly, we need to be careful with the IR region when Q' is sent to zero, since in this 
case the bulk to boundary propagator for the field dual to the outgoing on-shell photon 
spreads over the whole AdS space, therefore probing the IR region. In the particular case of 
DVCS, for which the outgoing photon is on-shell, we shall see in the next section that the 
Q' — > limit of the amplitude is well defined because we only need to consider transverse 
incoming and outgoing polarizations, so we do not need to be that careful in this respect. 
More generally, however, we must impose Neumann boundary conditions at the IR wall, 



which means that the Bessel functions defining (56) should be replaced according to 



K (Q'z) -> K (Q'z) - I (Q'z) , 

K^Q'z) -> KxiQ'z) + h(Q'z) . 



(63) 



The amplitude defined in this form is well defined, even when Q' becomes small. The same 
procedure should also be followed for the incoming off-shell photon, particularly when Q 
approaches the QCD scale (this will only happen for a reduced number of data points in 
section |5| . 

Thirdly, we also need to include the effects of the hard wall in the pomeron propagator 



[1], as briefly described in section 5.4 



3.2.1 Normalization 

We need to correctly normalize the function ty ab in the impact parameter representation for 



the hadronic tensor in (60), just as the impact parameter representation for the correlation 



functions is normalised according to (55), (56) and (58). Again, this can be done be equating 



the disconnected pieces of ( 59 ) and ( 60 ) . For ( 59 ) we have 



WSLfa) = (h\h)iC I d'ye*™ 



(64) 



where we note that (k^ka) = (2n) 3 y/s 5^(k 2 + k 4 ). To recover this expression from (60) 



one needs to follow the computation presented in an appendix of [T7]. This computation 
basically involves the integral representation of the bulk to boundary propagator derived 
directly from the CFT Regge theory, without making any reference to AdS space. The 



final result is exactly as in (55) and (56) for the four point correlation function, with the 



normalization of the current two point function given in (58). For details we refer the reader 
to reference [T7] . 
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It remains to fix the normalization constant C . The electromagnetic current operator 
has the well known form 

j a = e VfPflaiJf , (65) 
/ 

where the sum is over quark flavours, e is the electron charge and qf the corresponding quark 
charge in units of e. We shall be working in the kinematical region (47) already discussed in 
section 2.4 At weak coupling we can compute the two point function of this current, taking 
into account only the quark flavours that couple the photon probe to the target. Following 
the discussion in section [274 , in the kinematical range considered in this paper only the u, 
d, s and c quarks contribute. With this assumptions, the constant C computed using free 
Wick contractions in QCD is 

S 1 fl 

C = e 2 -^ — = 7.845 x 10" 4 . (66) 
4vr 4 9 v ; 

The above is a weak coupling computation. When fitting to experimental data, we are 
actually in a region where the coupling is not small. More specifically, we will obtain for the 
value of the intercept jo ~ 1-2 — 1.3, which clearly is in a region of finite coupling. Moreover, 
both the conformal pomeron and hard wall pomeron models consider an expansion around 
strong coupling. Thus we do not expect the normalisation constant in the current two 
point function to be given by (66) above. It should be given by a number of the same 
order of magnitude, but the exact numerical factor is a strong coupling computation that 
could only be done if we knew the dual to QCD. We therefore leave the constant C as a 
fitting parameter in what follows. The resulting value for C is indeed of the correct order 
of magnitude, confirming our expectations. 



3.2.2 Polarization Basis 

To compute the cross section for DVCS, and also for timelike deeply virtual Compton 
scattering where an outgoing timelike photon is emitted and then converted to a lepton or 
a qq pair, we need to contract the hadronic tensor with the photon polarization vectors. 
Let n\ and n' x be, respectively, the incoming and outgoing polarization vectors, normalized 
such that n 2 = n' 2 = 1 for transverse polarization, and n 2 = n' 2 = ±1 for longitudinal 
polarization (time-like or space-like photon). We can use the gauge freedom to impose the 
conditions 

n ■ h = , n' ■ k 3 = . (67) 

Throughout the computation we work with spacelike off-shell photons and then take the 
on-shell limit, or analytically continue for the time- like case, if necessary. 
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In light-cone coordinates (+, — , _L), the polarization of the transverse incoming photon 

is 

n A = (0,0,e A ), (A = 1,2) (68) 

where e A is a orthogonal basis of unit vectors on M 2 . The incoming longitudinal photon has 
polarization 

n 3 =i(vS,^=,0) , (69) 



where we define Q = yQ*. Note that all the rij are orthogonal (and unit normalized). 
For a transverse outgoing photon (on-shell and off-shell) we have 

n A = (o,2^,e' A ) , (A = 1,2) (70) 

where e' A is a orthogonal basis of unit vectors on IR 2 . For a longitudinal outgoing photon we 
have 

1 / r Q' 2 + ll 



n '* = Q>V S '^^- (71) 



Again all are orthogonal (and unit normalized). 
3.2.3 Deeply virtual Compton scattering 

We are finally in position to compute the DVCS cross section. We need to square the ampli- 
tude, and then average over the incoming transverse photons and sum over final transverse 
on-shell photons. For both incoming and outgoing transverse photons (A, A' = 1, 2) a simple 
computation shows that 

W xy = (n x ) a W ab (n' x ,)l = (e x ■ e£) W l = 5 XX > W l , (72) 

where we define 

W n = 2isQQ' J dl ± J J ^ * n (z) <S>(z) [l - e^ s ^] . (n = 0, 1) (73) 

It is now clear why the limit Q' — > of the amplitude Wxx' is well defined and erases any 
IR dependence on the outgoing photon wave function, since when computing W\ one takes 

Q' (ki(Q'z) + h{Q'z)^ l/z , (74) 



where we are accounting for the hard wall modification explained in (63). Finally we can 
square each amplitude, take the average over the incoming polarizations and sum over the 
final ones. We obtain 

\J2\ W ^\ 2 = \ W ^- ( 75 ) 

We remark that, in the longitudinal-transverse case we obtain after some algebra that 
Wlx' — 0, so that helicity is conserved, as expected. 
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3.2.4 Timelike deeply virtual Compton scattering 

Another interesting possibility is when the outgoing photon is timelike. It may then decay 
to a lepton pair (ee, fift or rf), or to a qq pair, provided the photon momentum is above the 
corresponding threshold. In particular, the process for lepton pair production is very clean 
and entirely within QED perturbation theory. Although we are not analysing in detail these 
processes in this paper, let us write for completeness the corresponding matrix elements in- 
volving longitudinal outgoing photons. In this case Q' = a/ Q' 2 is purely imaginary, so we 
need to be careful in analytically continuing the above expressions. Again helicity is con- 



served, so in addition to ( 72 ) we have the following non- vanishing longitudinal- longitudinal 
matrix element 

W LL = (n 3 ) a W ab (n' 3 )l = -W , (76) 



where W is also defined in (73). It would be interesting to further explore this process 



using AdS/CFT and pomeron physics. 



4 7*p — > jp cross section 

Our goal is to predict the experimental data collected at HERA for the differential and total 
cross sections in DVCS. For the differential cross section we have 

dcx . „n > 1 l^„ ir l9 \Wi\ 2 

which depends on the three kinematical invariants s, t and Q 2 . For the available data the 
total cross section 

a(x,Q 2 ) = ^| dt\ Wl \ 2 , (78) 
is obtained from integrating the differential cross section from to —1 GeV 2 . 



Let us then further simplify the expression for the function W\ in (73). First we do the 
angular integral in the impact parameter space 

r z ° d7 H7 f°° 
W n « 4vm QQ' — — * n {z) / dl I J (ql) [l - e 1 ^} , (79) 
Jo z z Jo 



where q = \q±\ = y—t and I — \l±\. Since we use a delta function for <&(z) we can further 
simplify this expression to 

f z " dz f°° 

W n « 47ris QQ' / -3 * n (z) / dl I J (ql) [1 - e^ s ' L) ] , (80) 
Jo z Jo 



where L = L(z, z*,l) defined in (52) is now computed at z = z*. 
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We will now consider three models for the phase shift function x{S,L): the AdS black 
disk, which is valid in a limited kinematical region of available data but very simple; the 
conformal pomeron at strong coupling that neglects confinement corrections to the reggeon 
trajectory but is still rather simple and reproduces data in a very satisfactory way over a 
large kinematical range; the hard wall pomeron that includes the effects of confinement in 
the previous model. 



4.1 AdS Black disk 

We shall consider a black disk model defined by a phase shift in the impact parameter 



representation of the hadronic tensor (60) given by 

1 - e ix{s ' L) = e(L s (S) - L) , 

where the saturation radius L s of the disk increases with energy as 

L a (S)*ulogS. (81) 

This model reproduces quite successfully DIS data at low x, in the kinematical region for 
which lmx(S, L) ^ 1 and for a probe above the QCD scale, i.e. for Q ^> Kqcd 13 EES]- In 
this simple model the phase shift is characterized by the single parameter u > 0. As long 
as the imaginary part of the phase shift grows as L decreases, one expects this model to be 
a good approximation provided the external kinematics imply that the dominant region of 
integration in the holographic variables z and z is such that L < L s . This model can be 
motivated using the BFKL pomeron at weak coupling, the gravi-reggeon and the phenomena 
of geometric scaling observed in DIS at low x, as explained in [19] . 

For the black disk model, the integral over the impact parameter can be better done by 
moving to the AdS impact parameter L where the phase shift is trivial. However one needs 
to be careful with the limits of integration. For fixed z, the upper limit of integration in 
L is determined by the size of the black disk L = L s (zz*s). The lower limit is determined 
by setting l± = 0, i.e. for L = \hx(z*/z)\. This region of integration, as a function of z, is 
shown in figure |1J However when changing integration variables (z, I) to (z, L) the region 
of integration in z is also bounded by the condition L s > |ln(^/ ^r) | , which introduces a 
lower bound of integration in z. The lower bound Z- is simple to understand from the AdS 
perspective. For a fixed target at z*, as the probe gets close to the boundary, the AdS 
impact parameter grows until it reaches the boundary of the disk for L(z, z*, 0) = L s (zz*s). 
Thus the condition defining z_ is given by 

— + — = (z„z*sY + (z_z*s)-" . (82) 
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Figure 4: The upper and lower bounds of integration in the H3 impact parameter L. Since the 
probe wave function is localized near the boundary, the dominant region of integration is within 
the interval z_ < z < z*. 



There are two solution of this equation, obtained one from another by u — > —uj. We are 
interested in the solutiorE] 

(z-z.s) u = — => zl + " = z^s)-" . (83) 
This lower bound is also shown in figure |4| We conclude that 

W n « AirisQQ'z* / ~^^ n {z) / dLsinhL J (ql) , (84) 

where the function 

/ = l(L, z, z*) = \J Izz* cosh L — z 2 — zl (85) 



is trivially obtained from (52). The L integral can be done by writing u = (ql) , with the 
result 

s f z ° dz 

W n « Am -QQ' — V n (z) l s J x {ql a ) , (86) 
Q J z_ z 

where l s = l(L s ,z,z*). Taking the Q' — > limit, and replacing for the expression of \&i we 
finally obtain 

27T 3 s r z ° 

W^-C — i-Q dzKr^lsJ^qls). (87) 
3 Q J z- 

This expression depends on the three parameters C, u and z*. In general it also depends 
on z , but for a probe localized in the UV region, this dependence drops out. In the next 



5 The other solution gives an upper bound for the z integration, corresponding to the probe placed deep 
in the infrared with respect to the target. We do not need to take this into account because the cut-off Zo 
will be close enough to z* . Also, since the probe wave function will be localized near the boundary in the 
UV, the result will in practice be independent of zq. 
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section we will use this expression in the differential cross section (77) and in the total cross 
section (78) to compare with HERA data. 



4.2 Conformal pomeron 

The next model we will consider for x{S, L) is a model for pomeron exchange in a conformal 
theory. In the case of an amplitude dominated by single pomeron exchange, the phase shift 



X = x(S, L) is small and is related to the conformal amplitude introduced in section 2.1 by 



B = 1 — e %x ~ —%X- Explicitly our single conformal pomeron exchange model is obtained 



by using expression (73) for W\ with the conformal amplitude computed in (41). 

Our conformal pomeron exchange model depends on three parameters. One of them is 
the position of the target z*, and is the same as in the previously discussed black disk model. 
The next parameter is p, related to the pomeron intercept via p = 2 — jo- m terms of the 't 
Hooft coupling, defined by the ratio i? 4 / a' 2 , it is given by p = 2/y/x. Finally, the amplitude 
will also depend on the overall constant Cg 2 , where C is given by the normalization of the 
current two point function and g\ by the impact factors of the external states. It is clear 
that for the single pomeron exchange we can not determine both constants independently. 

We expect the conformal pomeron model to give good results in the limit of Q 2 large, 
with x small, but not too small. We need x to be small in order to be in the kinematical 
regime where the pomeron exchange dominates. However, taking the limit of x — > of 



equation (41), we can see that its asymptotic behavior for the cross section is 



a~(l/x) 2 ^'- 1 )»A Q 2 7D ln 2 (l/x), (88) 
which violates the Froissart bound. This limit would also make Im% 3> 1, which means that 



at very small x (41) is no longer valid, and we expect multi-pomeron exchange effects to 
play an important part. If we perform the eikonal resummation of pomeron exchanges our 
amplitude then becomes the black disk model of the previous section, as shown by figure 



3j of section 2.3 However, the black disk amplitude (87) also violates the Froissart bound. 
This is not surprising, and it is not a contradiction, because both models are constructed 
from a conformal amplitude B(S,L), without a scale to define the Froissart bound in the 
first place. 

To better understand the amplitude as we go to very small x at fixed Q 2 , we need to 
realise that confinement effects also start to play a crucial role, as shown in [21] , and further 



discussed in section 4.3 below. Moreover, when Q 2 or the momentum exchange t approach 
the QCD scale, we also expect confinement effects to become important. This brings us to 
the hard wall pomeron model that we now discuss. 
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4.3 Hard wall pomeron 

The hard-wall model is obtained by placing a sharp cut-off on the radial AdS coordinate 
at z = zq. This changes the boundary conditions of the differential equations, and hence 
leads to novel results. The advantage of this model is that it is relatively straightforward 
to implement in AdS/CFT, leading to equations that can often still be solved analytically 
The cut-off also sets a scale in the theory, and therefore breaks conformal invariance and 
leads to a theory of confinement. Far away from the cut-off, its effects should be only 
weakly felt though, and we would recover once again approximate conformality, as in QCD. 
To see intuitively why the hard-wall model is confining, we can start by thinking about 
how the extra coordinate arises in the AdS/CFT correspondence in the first place. We 
can think of the graviton on the AdS side as being a product of two gluons from the CFT 
side. The radial coordinate in such a picture corresponds to the separation between the two 
gluons [50]. Therefore putting a maximum distance on this coordinate gives us a maximum 
separation between the gluons, and hence confinement. 

Let us commence by looking more closely at the relationship between the pomeron in the 
conformal and the hard-wall models. As mentioned before, the pomeron propagator in both 
these models satisfies the same differential equation, but with different boundary conditions 



jl]. First notice that at t — 0, equation (41) can be integrated in impact parameter l±, with 
the result 



POO 

x(r, i = 0, z, z) = Ikzz J dL sinh L x(r, L) 

J\\n(z/z)\ 



<\Mz/z)\ 

(M*/z)) 2 

= Z7T gl (cot (^) + i) (zz) "-^-^ , (89) 

where we defined the rapidity r = ln(a'S'). Similarly, the t — result for the hard- wall 
model can also be written explicitly 

Xhw(r,t = 0,z,z) = x(r,0,z,z) + J"(r, z, z) x(t, 0, z, z$/z) . (90) 

We see that at t — the hard-wall result can be represented as a sum of two conformal 
amplitudes, with the second one representing a conformal kernel on the other side of the 
hard-wall. The function 

^(r,z,^) = l-4 v ^ e " 2 erfc(r ? ), rj = ~ lo g(^/j)+4r 

V4r 



is set by the boundary conditions at the wall and represents the relative importance of the 
two terms and therefore confinement. This function varies between —1 and 1, approaching 
— 1 at either large z, which roughly corresponds to small Q 2 , or at large r corresponding 
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to small x. It is therefore in these regions that confinement is important. This relationship 
has been explored before in [24] . Since we are working at small x, the size of J 7 will roughly 
vary between —0.1 and —0.4 for the data here analysed. 



To use the hard- wall propagator in equation (80) we need to have the propagator at 
finite impact parameter I. Work is in progress in evaluating it numerically. In this paper, 
we used the same approximation as in [21]. It can be shown that at large / the eikonal for 
the hard-wall model has a cut-off 

Xhw(i~, I, z, z) ~ exp[-miZ - (m - mi) 2 I 2 /Apr] , (92) 

where mi and tjiq are solutions of 

d z (z 2 J ( mi z))\ z=zo = 0, d z (z 2 J 2 (m z))\ z=zo = 0. (93) 

At small I we assume x still has the t = form of a sum of two conformal kernels, but now 
with I dependence as well 

Xhl(T,l,z,z) ~ Xc(t,1,z,z) + J r (T,z,z)xc(T,l,z,z 2 /z). (94) 

Finally we will introduce a normalization function, C(r, z, z'), independent of I, which en- 



sures that the result (90) is reproduced after integrating over /. Putting it together, the 



approximation we use for x 111 the hard-wall model is 

Xhw{r, I, z, z) = C(r, z, z) D(t, I) X^(r, I, z, z), (95) 

where 

n / n • A exp[-miZ - (m - m^P/ipr] \ 

D{t, Z) = mm 1, = „ — r (96) 

V exp[-m!2;o - (m - mi) 2 ^/4prjy 

is the exponential cutoff at large /. We expect this to be a good approximation in the 

range of (x, t, Q 2 ) where experimental data is currently available. This model has one more 

parameter, zq, than the conformal pomeron modelj^] 



4.4 Eikonalization 

As mentioned in the previous two sections, at very small x single pomeron exchange violates 
the Froissart bound, which is common to both the conformal and the hard-wall models. One 

6 mi 1.6/zo and m sa 3.8/zo are the masses of the lightest spin 2 glueballs on the lowest two Regge 
trajectories. For more on calculating glueball masses see [47] . 

7 Note that even though z$ appears in the equation for the conformal model as well, it is not really a 
parameter there since, due to the shape of the probe wave function, as long as Zq J> 5 GeV -1 the result will 
be independent of its value. 
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way to avoid this difficulty is to fit directly equation (80) to data, without expanding the 
exponent in e* x , where \ is taken as the amplitude for single pomeron exchange. This 
would introduce an extra parameter in the fit, the overall normalization factor. Staying at 
single pomeron exchange level, only the product of this parameter with g$ appears, and so 
it cannot be fit. However, in the kinematical region where data is available, the size of x 
is small, and eikonalization does not change the value of our parameters, nor improve the 
quality of the fit. This is consistent with the conclusions of |24j. As data at smaller values 
of x, with Q 2 fixed, becomes available, we expect eikonalization to play a more important 
part, and to enter the black disk regime. 



That being said, we will briefly revisit eikonalization in section |5.3[ Unlike for single 
pomeron exchange, C and g\ can be fit separately in the eikonal model, since g\ is exponen- 
tiated and C is not. This will therefore allow us to estimate the size of C and to compare 
it with that obtained in the black disk and with the weak coupling result (66). 

It is also worth noting that, in contrast with the conformal models, the hard- wall x given 
in equation (95) will satisfy the Froissart bound, due to the exponential cut-off, (96). 



5 Data analysis 

We will now compare our results to data collected at HERA. We will use the data from 
both ZEUS and HI experiments, taken from their latest publications on DVCS 



The data sets have results for both the differential cross section (77), and the cross section 



(78). Figure [5] shows all the data points analysed in this paper, as a function of S ~ Q/x 
and Q. Compared to DIS, the data sets are smaller, with 52 and 44 points respectively. In 
a subsequent publication, we plan to extend our study to vector meson production as well, 
which will give us more data. 

In turn, we will look at each of our models. We will see that the fits give us precisely 
what we expect from the theory. In the case of the differential cross section, the best fit is 
obtained using the hard-wall model. The conformal pomeron model begins to show signs 
of breaking down at small values of t, which we expect to be the case due to the effects 
of confinement, which is precisely what the hard-wall model incorporates. The black disk 
model gives a very good agreement with experiment in a limited regime where we a priori 
expect it to work - when the size of x starts increasing, but before we reach the region where 
confinement is important. 

For the cross section, the values of x and Q 2 are in the region where we expect the 
conformal model to work, and before confinement starts becoming important. This is borne 
out in the fits, with both the conformal and hard wall models giving very good fits. The 
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10 10 10 10 

Q Q 

Figure 5: (Left) Data points to fit the differential cross section. Note that for each point we have 
4 different values of —t, so we have a total of 52 points. In the case of the black disk model we 
restrict the analysis to the filled dots (20 points), as explained in the text. (Right) Data points to 
fit the total cross section, giving a total of 44 points. The filled dots represent the points we select 
to fit the cross section using the black disk model (10 points). In both figures Q is in GeV. 



conformal fit is slightly better, but the difference is not large. Once again, the black disk 
model also gives a very good fit in the region where we expect it to work. 

Note that the normalization of the photon wavefunction C is a parameter in all models, 
but for single pomeron exchange it only appears in the combination Cg\ and hence these 
two parameters cannot be fit separately. The results we present below in sections |5.2| and 



5.4 were obtained by fixing C to the reference value (66), given by the weak coupling 
computation. Therefore bear in mind that even though we quote a value for it should be 
understood that we only know the value of the product Gg\. 



5.1 Black disk 

The black disk model is expected to be valid in a kinematical window where QCD is ap- 
proximately conformal and where the amplitude for single pomeron exchange becomes large. 
This reduces the number of available data points to a small subset. The criteria to select 
data points was discussed in detail in the original work [TJ, here we only summarise it. 
There are three conditions we need to impose: {%) Q > Q m in ~ 2 GeV, so that our probe 
does not feel confinement effects; (ii) To be inside the saturation region, defined by the 
condition Imx ~ 1, we need to impose uln(z*Q/x) > ln(z*Q), which can be derived from 
the expression of the conformal amplitude for the exchange of a reggeon; (Hi) We need to 
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be in the Regge limit of large s which, given the value of the other kinematic invariants, 
gives the condition z*Q/x > 10 3 . The straight lines plotted in figure [H] correspond to these 
three conditions, so that we are left with 20 points for the differential cross section and 10 
points for the cross section. 

The parameters we get from the differential cross section data are 



u = 0.243 ± 0.045 , z» = 3.50 ± 0.89 GeV" 1 , C = 0.0013 ± 0.0010 , (97) 
which gives us a x 2 of 

XL./. = 0-97, (98) 

showing that the black disk fits in the range we expect. Analogously, the parameters from 
the total cross section data are 

w = 0.275 ±1.581, z m = 1.39 ± 124.80 GeV -1 , C = 0.00094 ± 0.00530 , (99) 

corresponding to a x 2 °f 

Xlo.f. = 1-22. (100) 

In figure [6] we see the plots for these fits. Considering that we are fitting a very small 
number of data points, the results are very good. Both values for the exponent u are very 
close, and correspond to an intercept j = 1 + u between 1.2 and 1.3, as expected. Also, 
for the overall constant C and the proton scale we obtain values with the expected order 
of magnitude, as given by the weak coupling result (66) and the proton mass, respectively. 
Although results are making sense, with very good values for x 2 , a word of caution is required 
for the very large errors obtained in the estimated parameters. This happens because we 
have a very small number of data points in a very limited region where the black disk 
model only starts to describe the physics. With data for smaller values of x, inside the 
region defined by the black disk criteria plotted in figure |5j we would expect to be able to 
considerably reduce these errors. 



5.2 Conformal pomeron 

The parameters we fit for the conformal pomeron are the coupling to external states g$ (in 
fact C^ofl t ne ra dial position of the proton z* and p, which is related to the intercept by 
jo = 2 — p. For the differential cross section the values we get are 

£o = l-95 ±0.85, z* =3.12 ±0.160 GeV^ 1 , p = 0.667 ± 0.048 . (101) 

8 Throughout the paper we fit using a weighted non-linear least-squares method, and the parameter errors 
represent the 95 percent confidence interval for the parameter estimates. 

9 Note that the definition of is different than that in [21] , so the values for this parameter are not 
directly comparable, but they can be simply related to each other. 
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W = 70 GeV 



W = 82 GeV 




Figure 6: The fit of the black disk model to the data. The first three figure are for the differential 
cross section and the last one for the cross section. In these figures W is the center of mass energy, 
W = \/s, and we use units of GeV. 



32 



corresponding to a x 2 °f 

Xlo.f. = 1-33. (102) 

We can see in figure [7] the comparison of our model to data. Even though —t is larger than 
Aq CD for all values where we have data, at the smallest values it is close, and there might 
be some lingering effects of confinement felt. If we exclude the lowest value of \t\ from each 
graph, we get a better fit, with a x 2 of 

Xlo.f. = 0-76. (103) 

For the cross section, we fit the same three parameters. The values we get are 

#0 = 8.79 ±4.17, z* = 6.43 ± 2.67 GeV" 1 , p = 0.816 ± 0.038 . (104) 

Here the integration over t is likely to reduce the weight of the region of low t to the data, 
so that confinement effects are less important, resulting in a better fit with a x 2 °f 

Xlo.f. = 1-00. (105) 

In figure [7] we see how the model compares with data. We see that we have a very good 
fit. Note that in the figure we skip over some values of Q 2 in order to avoid cluttering the 
graph. 



5.3 Eikonalization 



As mentioned in section AA fitting the eikonal model to data does not improve our fits. 
Recalling that g\ in (101) is actually Cg$, with C given by the reference value (66), the 



value for the parameters we obtain for the eikonal fit does not change much from that in 
(101). The only difference is that the eikonal fit gives a g 2 ^ 0.05 and a value of C that 



compensates for that change, keeping Cg 2 , approximately constant. Indeed, by running the 



fit for the eikonal model the values of p and are the same as those in equations (101) 
and (104). Now it turns out that once g\ ~ 0.5 or lower, the eikonal fit does not change 



much because we enter the regime where the phase shift in the exponential is small, and the 
amplitude dominated by single pomeron exchange. It is nevertheless encouraging that the 



estimated value for C has the order of magnitude of (66 ). The very same order of magnitude 



for the constant C was also obtained in the black disk fitting (97). In both cases there are 
large uncertainties in the estimated value of C. 

It is also useful to analyse the size of % f° r single pomeron exchange in the kinematical 
region of validity of the black disk model, figure [8] Due to the aforementioned uncertainty in 
C, it is not the exact values for x that are important, but rather the pattern that emerges. 
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W = 40 GeV 



W = 70 GeV 




Figure 7: The fit of the conformal pomeron model to the data. The first five figures are for the 
differential cross section data and the last one for the cross section. W is the center of mass energy 
and we use units of GeV. To avoid cluttering the last figure we did not plot all of the Q 2 values. 
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10" 10' 



Figure 8: Contour plot of the phase shift x- Q m GeV. 



We see that the lines of equal x are roughly parallel to the line \og(Qz*), as expected, and 
the black disk fits in the range where x grows, excluding the points with a small value of %■ 
Also excluded is the region where we believe confinement will start to play a role. 
The same conclusions are supported from the analysis of the total cross section. 

5.4 Hard wall pomeron 

As mentioned earlier, the hard wall pomeron adds an additional parameter, the position zq 
of the cut off in AdS. The parameters we obtain by fitting are 

#o = 2.46 ±0.70, z* = 3.35 ±0.41 GeV -1 , p = 0.712 ± 0.038 , z = 4.44 ±0.82 GeV -1 . 

(106) 

corresponding to a x 2 °f 

XL./. = 0.51. (107) 

The fit is clearly better than the conformal one, as can also be seen from figure [9] The main 
improvement that it offers is taking into account the effects of confinement, which are still 
felt at low \t\. For the cross section (also in figure [9]), the values for the parameters we get 
are 

#0 = 6.65 ±2.30, z* = 4.86 ± 2.87 GeV" 1 , p = 0.811 ± 0.036 , z = 8.14 ± 2.96 GeV -1 . 

(108) 

corresponding to a x 2 of 

xL.f. = 1-03. (109) 
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W = 40 GeV 



W = 70 GeV 




Figure 9: The fit of the hard- wall pomeron model to the data. The first five figures are for the 
differential cross section data and the last one for the cross section. W is the center of mass energy 
and we use units of GeV. To avoid cluttering the last figure we did not plot all of the Q 2 values. 
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Based on the fit to DIS data |24| . we expect the effects of confinement only to become 
important once Q 2 is about 1 GeV 2 or lower, and above those energies the conformal and 
hard-wall models to give very similar results. The above fits confirm these expectations, 
both with the x 2 values and the fact that the fitted paramaters for both models are close 
to each other. 



6 Conclusion 

One of the major goals of particle physics is to accurately calculate and predict cross sections 
for scattering experiments. Therefore novel techniques that allow us to do such calculations 
are always welcome, both from a practical point of view, and a theoretical one, since they 
deepen our understanding of particle interactions. Gauge/string duality already has many 
interesting applications, and in this paper we have carried out further steps in using it to 
study diffractive processes in QCD. 

Indeed, pomeron exchange is the dominant contribution to the cross section in low x 
physics. When all scales are above Aqcd, we can explore conformal Regge theory and 
compute the cross section arising from one reggeon exchange. This leads to the confor- 



mal impact parameter representation (51) for the transition amplitude T a b(ki), where the 



conformal amplitude B(S, L) can be computed both at weak coupling M3J) and at strong 



coupling (41). The approximate conformal symmetry in processes dominated by pomeron 
exchange, when the amplitude grows beyond the single pomeron exchange, then leads to the 
AdS black disk model proposed in [7j and successfully tested against DIS data inside the 
saturation region. To extend the analysis to a larger kinematical range one could consider 
either the weak or the strong coupling expansion for reggeon exchange, and then include 
confinement effects that break conformal invariance. In this paper we considerer the strong 
coupling expansion for the exchange of a pomeron, first identified by Brower, Polchinski, 
Strassler and Tan in |3] as the graviton Regge trajectory in AdS, or the BPST pomeron. 
Its advantage, as compared to the perturbative BFKL pomeron, is that it unifies the soft 
and the hard pomerons through diffusion in AdS space, which was also shown explicitly in a 
comparison to data in [24J. This allows for comparison in a much larger kinematical range, 
unaccessible to traditional perturbative QCD methods, including for energy scales close to 
Aqcd, where the processes are inherently non-perturbative. The pomeron's universality in 
a very wide range of high energy scattering processes then makes such a study very useful. 

After reviewing and further elucidating some of the properties of the conformal pomeron 
Regge theory, we computed the hadronic Compton tensor and considered three models: the 
AdS black disk model, the conformal pomeron model and the hard wall pomeron model, 
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which incorporates confinement. Some of the computations have already been done when 
comparing to DIS data for the black disk model [7] and both pomeron models [24J, but we 
made several improvements necessary to compute the amplitude when t ^ 0, as required 
for deeply virtual Compton scattering. This is the natural next step after DIS, due to the 
similarity between both processes. 

In section [5] we carried out the data comparison for the above three models. We found a 
very good agreement with data, where each of the models is applicable. We also found that 
the parameters are in the same range as those found in the previous works [7J 121] , with the 
parameters for the cross section being the closest. This is natural, since DIS also measures 
the total cross section. Furthermore, all the main previous conclusions are supported: inside 
saturation one observes a conformal (AdS) black disk j7]; the conformal pomeron fits the 
data well in the conformal region, the hard wall model becomes more accurate when Q and 
t approach the QCD scale in a region where confinement effects become important before 
saturation starts to occur ET 



Our analysis in sections and <xA further suggests that confinement effects become 
important when —t is of order 0.1 GeV 2 , and the hard- wall model can explain this data as 
well, further cementing its status as a useful tool in studying the string duals of confining 
theories, such as QCD. 

Furthermore, as for DIS, the fit for the pomeron intercept again has a value in the 
crossover region between strong and week coupling, jo = 1-2 — 1.3. The value we obtain 



for the intercept, as a function of the 't Hooft coupling is represented in figure 10, together 
with the corresponding weak and strong coupling curves. Since we are in an intermediate 
region, it is very possible that the weak coupling expansion is as good as the strong coupling 
one. This raises the question whether one can construct a weak coupling hard wall model. 
Indeed, an attempt to fit DIS data using the BFKL pomeron with a IR cut-off appeared in 
[ST] , but the strong coupling model fits data better in the IR region. It is therefore possible 
that using the (stringy) AdS dual as a tool to introduce a IR cut-off, just as in the strong 
coupling case, will give an equally good fit to data. We leave this question to future research. 

Our goal is to develop a framework for a study of a wide range of high energy scattering 
processes, and eventually use all these processes to constrain the parameters of our model. 
For example, other processes where we can apply an analogous study are vector meson 
production and also proton-proton total cross sections. It is also important to move beyond 
2—7-2 scattering. Indeed, work on applying pomeron exchange in AdS to double diffractive 
Higgs production, a 2 — > 3 process, should appear shortly [5^] . 

We have taken here an important step for such a unified framework, but though much is 
taken, much abides. As a next step we plan to extend this study to vector meson production, 
where the coupling of the mesons to the pomeron is determined by the IR physics deep 



38 



Jo 



2 



0.5 



1.5 



1 




2 



4 



6 



8 



aN 



Figure 10: The comparison of the pomeron intercept for weak and strong coupling. The green and 
blue dotted lines give the dependence of the pomeron intercept to first and second order on the 
't Hooft coupling aN , respectively. The solid line gives the strong coupling dependence obtained 
from strings in AdS. Figure taken from reference [3J. 

inside AdS. This will be an important study, among other reasons for allowing a better 
examination of the interplay between confinement and saturation. Another useful project 
is to numerically calculate the hard-wall pomeron propagator in impact parameter space, 
thus improving the approximation we used here. It would likewise be quite interesting to 
use a different model of confinement, allowing us to identify which of its features are model 
independent. We expect to return to such questions in the near future. 
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A Graviton exchange in AdS 



The impact parameter representation (51) assumes conformal symmetry and Regge kine- 



matics. In this section we compute this amplitude using the AdS / CFT duality, in particular 



we check the non-trivial form of the tensor function (55) associated to the coupling of the 
vector current (dual to a U(l) AdS gauge field) to the pomeron (dual to the graviton Regge 
trajectory). To fix our notation we start by writing the action for the gravity fields 

S = 1 t 2 I d " x ^~ 9 ( n ~ 12 ) + Sm ■ ( U °) 

In our units the AdS vacuum has radius 1. In the context of M = 4 SYM the gravitational 
coupling k = 2tt/N and the AdS radius R A = 4irg s Na' 2 . For the matter fields we shall 
consider a U(l) gauge field A, dual to the electromagnetic current operator, and a complex 
scalar field with mass given by m? = A(A — 4), with a bound state dual to the target 
hadron (after introducing a gauge/gravity model for confinement). The matter fields action 
is therefore given by 

S m = J d 5 xV^ (-\ F 2 - i(«90) 2 - ^A(A - 4)/) , (111) 

where F = dA. 



A.l Propagators 

We denote the bulk to boundary propagator for the scalar field by IT(X, y) where X = (z, x) 
are the Poincare coordinates of the bulk point and y the coordinate of the boundary point, 
with x, y G M 4 . We have 

with w = x — y. It will be useful below to have the Fourier transform of this propagator 

n(z, k) = j d A w e* k - w IT(X, y) = f |^y Q A ~ W 2 _ A (Q^) , (113) 
where Q = y/k?. Taking the boundary limit we have z A-4 II(X, y) — > 1 • 8{x — y). This 



normalisation fixes the constant C in (58) to 



- A -2 T(A + 1) , . 

C = 2 ^^K^2 Y (114 > 
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For the bulk to boundary propagator of the U(l) gauge field Il Aia (X, y), working in the 
Lorentz gauge D^Il^a = 0, we have 

n fY \ 1 (3^ 2 - w 2 ) r] ab + 4w a w b 

Tl ba (X,y) = — — - — — , U za {X,y) = 0, (115) 

2ir z (z z + w z ) a 

where again w = x — y. The Fourier transform is then 

Tl b a(z,k) = J d A we ik - w U ba (X,y) = (r) ab - (Qz^Qz) . (116) 

It is now trivial to see that taking the limit near the boundary the propagator generalises 
U ba (X, y) — > rjab • S(x — y), such that the Lorentz gauge condition d b U ba = is satisfied on 



the boundary. Computing the two point function we obtain a normalisation in (58) with 



c=^- ( 117 ) 



Let us finally consider the bulk to bulk propagator for the graviton. In the Regge limit 
the propagator equation for a spin 2 and dimension 4 field simplifies to (see [HJ [53] for 
details) 

(n + 2)n a ^„(X,X) = ig a{tl g v) p5 AdS (X,X) + ■ ■ ■ . (118) 

where • • • represent terms that are not important in the Regge limit. As usual, given the 
energy-momentum source T^, the metric perturbation is given by 

(119) 



h a0 (X)=in J d 5 X^U a ^ u (X,X)T^(X) 



In the following discussion we will not need the full form of the graviton propagator. Instead 
we will only need to know the integrated propagator 



J dx + dx-U ++ __(X,X) = ~2i(zz)- l Yi ± (L) , 



(120) 



where we use light cone coordinates x = (x + ,x ,x±) on the M 4 boundary and similarly for 
x. n_i_(L) is the scalar propagator in H% with mass squared 3, and it depends only on the 



invariant L given in (52) with l± = x± — x±. In other words, Il_|_ satisfies the Euclidean 
equation 

(n HB -3)H ± (x,x) = -6 Ha (x,x), (121) 
where x = (z,x±) and x = (z,x±) in Poincare coordinates on H%. Explicitly we have 



1 e~ 3L 
2tt Y^e- 



^{L) = -^—^. (122) 
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A. 2 Witten diagram 

We have now the necessary ingredients to evaluate the Witten diagram that computes the 
correlation function (50) in the Regge limit, which is dominated by the t-channel graviton 
exchange. Feynman rules in AdS give for G ab (yi) = (j a (y i) O (2/2) jb(y '3)0(1/4)) the expression 

'd b X d 5 X 



A a b(yi 



[IK] 



z 5 z 5 

„4-4 



2d h U a]a (X, yi ) g a ^(z) 2d h U fi]b (X, y 3 ) 



(16zV) n ++ __ (X, X) «9+IT(X, y 2 ) d + Tl(X } y 4 ) 



(123) 



where we used that the graviton couples to the energy momentum tensor T^ v as given in 
(119), and we note that the factor lQz 4 z 4 appears from metric elements g + ~ computed at z 
and z. Evaluating the Fourier transform 

A a b(ki / 



I (n dv - **) A ab ( Vi ) 



(124) 



we obtain 

Aab (ki 



d 5 X d^X 

r e-^ +k *>* 2k l]r U a]a (z, k x ) g^(z) 2* S[ _n fl6 (z, k 3 ) 



z 5 z'° 



(16z 4 ^ 4 ) n ++ __(X, X) e -^ k2+k ^ s k 2+ Tl(z, k 2 ) h + U(z, k, 



4) , 



(125) 



where the momentum space bulk to boundary propagators are given in (113) and (116). 
Next we insert the external kinematics, as given in (48) and (49). In the Regge limit the 
result simplifies to 

dz dz 

' ' - 2k n _U a]a (z, h) g a \z) 2k 3[ _ U m (z, k 3 ) 



z 3 z 3 



di3y do^ 



(I6z 2 z 2 ) n ++ __(x,x) 



(126) 



A ab (k i )=-K 2 (2<K) 4 5(0) J dl x e^ 4 
k 2+ Yi(z, k 2 ) k i+ Tl(z, k A ) 

where l± = x± — x± . We replaced the space-time volume by the momentum conserving 
S function, V = (27r 4 )<5(0), to extract the amplitude T ab introduced in (50). The integral 
inside the large parenthesis gives precisely the scalar propagator in H 3 given in (120). Thus, 
we obtain finally 

dz dz 



Tab(ki) 



dh e iq±4A 



z 3 z 3 



z" (t ab V + uM - Tl(z, k 2 ) U(z, h)(16zz) U ± (L) 



(127) 



where is defined in (1561) and 



tab = (Q 2 V~a - h-kia) (Q' 2 V-b - k 3 -k 3b ) , 

Uab = QQ' (h-k 3 ^r] ab - k X -k 3a r}„ b - k lb k 3 ^ a + (k x ■ k 3 ) r/_ Q ?7_ 6 ) . (128) 
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A simple computation shows that 

/ s Q' 2 + ql v^?i \ 



t 



ab 



Q 2 \Q 2 {Q' 2 + ql) jsQM 



V 







ab 



/ 



( \ 


\0 f s q\ h) 



(129) 



Thus we recover the tensor function defined in (55) 

4z 4 



(130) 



with C = 6/n 2 . Also, using C as given in (114) we recover the scalar function in (57), more 
precisely 

<^>{z) = U{z,k 2 )U{z,k 4 ). (131) 



Finally we can read the strong coupling phase shift associated to graviton exchange 



X (S,L) = - Y SU ± (L), 



(132) 



in agreement with j44J [45] . This computation is a non-trivial check of the impact parameter 
representation introduced in [39], which only relies on conformal symmetry. 
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